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Abstract 

Joint power and admission control (JPAC) aims at supporting a maximum number of links at their 
specified signal to interference plus noise ratio (SINR) targets while using a minimum total transmission 
power in the interference-limited network. Almost all of the existing works on JPAC assume perfect 
channel state information (CSI). However, this kind of assumption generally does not hold true. In 
this work, we consider the chance (probabilistic) SINR constrained JPAC problem, where a random 
channel model is used and each link's SINR outage probability is enforced to be not greater than 

> 

a specified tolerance. The chance SINR constraints are computationally intractable in general even 
though the channel distribution information (CDI) is assumed. To circumvent the extreme difficulty 
of chance SINR constraints, we use sample (scenario) approximation scheme to simplify them. The 

(N 

sample approximation scheme not only transforms the difficult chance SINR constraints into finitely 
many (dependmg on the sample size) simple Imear constramts, but also relaxes the assumption of CDI, 
^ thus being distributionally robust. 

Furthermore, we reformulate the sample approximation of the chance SINR constrained JPAC prob- 
ed lem as a group sparse minimization problem and then relax it to a second-order cone program (SOCP). 

The solution of the SOCP relaxation problem can be used to check the simultaneous supportability of 
all links in the network and to guide an iterative link removal procedure (deflation). Instead of solving 
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the SOCP with general-purpose solvers like CVX, we exploit its special structure and custom-design an 
efficient algorithm. We illustrate the effectiveness and efficiency of the proposed sample approximation- 
based deflation approaches by using a nonrobust counterpart as the benchmark. 

Index Terms 

Chance SINR constraint, distributionally robust, group sparse, power and admission control, sample 
approximation. 

I. Introduction 

Joint power and admission control (JPAC) has been recognized as an effective tool for 
interference management in cellular, ad hoc, and cognitive underlay wireless networks for two 
decades Pl-p9|- The goal of JPAC is to support a maximum number of links at their specified 
signal to interference plus noise ratio (SINR) targets while using a minimum total transmission 
power when all links in the interference-limited network can not be simultaneously supported. 
JPAC can not only determine which interfering links must be turned off and rescheduled along 
orthogonal resource dimensions (such as time, space, or frequency slots), but also alleviate the 
difficulties of the convergence of stand-alone power control algorithms. For example, a long- 
standing issue associated with the Foschini-Miljanic algorithm [3| is that, it does not converge 
when the preselected SINR levels are infeasible. In this case, a JPAC approach must be adopted 
to determine which links to be removed. 

The JPAC problem can be solved to global optimality by checking the simultaneous support- 
ability of every subset of links. However, the computational complexity of this enumeration 
approach grows exponentially with the total number of links. Theoretically, the problem is 
known to be NP-hard to solve (to global optimality) and to approximate (to constant ratio 
global optimality) [[T|, Q, pO| , so various heuristics algorithms p|-p9| have been proposed. 
In particular, the reference [[T| proposed a convex approximation-based algorithm, called linear 
programming deflation (LPD) algorithm. Instead of solving the original NP-hard problem directly, 
the LPD algorithm solves an appropriate LP approximation of the original problem at each 
iteration and use its solution to guide the removal of interfering links. The removal procedure 
is terminated until all the remaining links in the network are simultaneously supportable. The 



reference 1 10| developed another LP relaxation-based new linear programming deflation (NLPD) 
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algorithm for the JPAC problem. In pO|, the JPAC problem is first equivalently reformulated as 
a sparse ^o-minimization problem and then its £i -convex relaxation is used to derive a LP, which 
is different from the one in [1|. Again, the solution to the derived LP can guide an iterative 
link removal procedure, and the removal procedure is terminated if all the remaining links in 



the network are simultaneously supportable. Similar ideas were also used in [13|, [19| to solve 
the joint beamforming and admission control problem for the cellular downlink network. 

The joint power/beamforming and admission control problem in the above-mentioned refer- 
ences except [[T| assumes perfect channel state information (CSI). However, such an assumption 
generally does not hold in practice [20|. In [[1], the authors considered the worst-case robust 



JPAC problem with bounded channel estimation errors. The key in [[TJ is that the relaxed LP 
with bounded uncertainty can be equivalently rewritten as a second-order cone program (SOCP). 
The overall approximation algorithm remains similar to LPD for the case of perfect CSI, except 
that the SOCP formulation is used to carry out power control and its solution is used to check 
whether links are simultaneously supportable in the worst case. 

In this paper, we consider the chance (probabilistic or outage-based) SINR constrained JPAC 
problem, where a random channel model is used and each link's SINR outage probability 
must be kept below a given tolerance. Unfortunately, this kind of optimization problems turn 
out to be extremely difficult to solve in an exact fashion, since chance SINR constraints are 
computationally intractable even though channel distribution information (CDI) is assumed. 



Thanks to recent advances in chance constrained optimization [21 1, |22|, the chance SINR 
constraint can be approximated via sampling with provable performance guarantee. The sample 
(scenario) approximation typically does not require explicit CDI. Instead, it only requires that 
the random channel parameters can be sampled. Therefore, the sample approximation technique 



is distributionally robust. As remarked in [23 1, most wireless systems are constantly sampling 
channel fading processes as part of the operation for coherent detection and adaptive modulation. 
Hence, the above sampling requirement can be easily satisfied. 

We further show that the sample approximation of the chance SINR constrained JPAC problem 
can be equivalently reformulated as a group sparse ^o-minimization problem, which is a nontrivial 



extension of the sparse formulation in [ 10 1. Then, we consider its -relaxation, which is an SOCP 
(different from the one in [[T|), and use the solution of the SOCP relaxation problem to check 
the simultaneous supportability of all links in the network and guide an iterative hnk removal 
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procedure (deflation). Instead of relying on standard SOCP solvers to solve the derived SOCP, we 
exploit its special structure and custom-design an efficient algorithm. Moreover, we develop an 
easily checkable necessary condition for all links in the network to be simultaneously supported at 
their desired SINR requirements. This condition allows us to iteratively remove strong interfering 
hnks and therefore significantly accelerate the deflation process. Numerical results show that the 
proposed sample approximation-based deflation approaches perform very well. 

Notations: Denote the index set {1,2,-- - ,K} by /C. Lowercase boldface and uppercase 
boldface are used for vectors and matrices, respectively. For a given vector x, the notations 
max{x}, min {x} , (x)^, and ||x||o stand for its maximum entry, its minimum entry, its /c-th entry, 
and the indicator function of x (i.e., ||x||o = if x = and ||x||o = 1 otherwise), respectively. 
The expression xi 0x2 represents the Hadamard product and max{xi,X2} (min{xi,X2}) the 
component- wise maximum (minimum) of two vectors xi and X2. For any subset X C /C, Aj 
stands for the matrix formed by the rows of A indexed by X. The spectral radius of a matrix 
A is denoted by p(A). We use (Ai, A2) to denote the matrix formed by stacking matrices Ai 
and A2 by column and use (Ai; A2) to denote the matrix formed by stacking Ai and A2 by 
row. Similar notations apply to stacking of vectors and scalars. Finally, we use e to represent 
the vector of an appropriate size with all components being one, I the identity matrix of an 
appropriate size, and the matrix of an appropriate size with all entries being zero except its 
A;-th column entries being one, respectively. 

IL Review of the NLPD algorithm 
The algorithms developed for the chance SINR constrained JPAC problem in this work are 



based on the NLPD algorithm [ 10 1 for the JPAC problem that assumes perfect CSI. To streamline 
the presentation, we briefly review the NLPD algorithm in this section. The basic idea of the 
NLPD algorithm is to update the power and check whether all links can be supported or not. 
If the answer is yes, then terminate the algorithm; else drop one link from the network and 
update the power again. The above process is repeated until all the remaining links can be 
simultaneously supported. 

Specifically, consider a J^-link (a link corresponds to a transmitter-receiver pair) interference 
channel with channel gains g^ j > (from transmitter j to receiver k), noise power T]k > 0, 
SINR target 7^ > 0, and power budget pk > for k,j E }C. Denote the power allocation vector 
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by p = {pi,P2, - ■ ■ ,Pk)^ and the power budget vector by p = {pi-,P2r " iPkY ■ Treating 
interference as noise, we can write the SBSFR at the fc-th receiver as 

SINR,(p) = , Vfce/C. (1) 

Vk + 2_^ gkjPj 

Correspondingly, we introduce an equivalent normalized channel. In particular, we use 

q= (9i,?2,--- (2) 

with qk = Pk/Pk to denote the normalized power allocation vector, and use c = (ci, C2, ■ ■ ■ , ck)'^ 
with Cfc = i'jkVk) /igk,kPk) > to denote the normalized noise vector. We denote the normalized 
channel matrix by A G M^^^ with the (/c, j)-th entry 

r 1, iffc=j; 

= \ lk9k,jPj 
I gk,kPk 

With these notations, it is simple to check that 



SINRfc(p) > 7fc ^ (Aq - c);, > 0. 

Based on the Balancing Lemma [[241, we reformulate the joint power and admission control 
problem as a sparse optimization problem 

inin ^^g^||(c- Aq)fc||o + ap^q 
*^ (3) 
S t. < q < e, 

where < a < ai := 1/e^p. In the above, e is the all-one vector of length K. Since problem 
(|3]) is NP-hard, we further consider its £i -convex relaxation (which is equivalent to an LP) 

min ||c — Aq||i + ap'^q 
^ (4) 

s.t. < q < e. 

By solving Q, we know whether all links in the network can be simultaneously supported or 
not. If not, we drop one link (mathematically, delete the corresponding row and column of A 
and the corresponding entry of p and c) from the network according to some removal strategy, 
and solve a reduced problem Q until all the remaining links are supported. 
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III. Problem Formulation 

Consider the chance SINR constrained JPAC problem, where the channel gains {gk,j} and the 
noise power {r]k} in the SINR expression ([1]) are now random variable^ Since {gkj} and {rjk} 
in ([T]) are random variables, we need to redefine the concept that a link is supported. We call 
link k is supported if its outage probability is below a specified tolerance e G (0, 1), i.e., 

P(SINRfe(p) >7fc) > 1-e, (5) 

where the probability is taken with respect to the random variables {gkj} and {r]^} . 

The chance SINR constrained JPAC problem can be formulated as a two-stage optimization 
problem. Specifically, the first stage maximizes the number of admitted links: 

max liSI 

P,S 

s.t. P (SINRfc(p) > 7fc) > 1 - e, A; G 5 C /C, (6) 
< p < p. 

We use So to denote the optimal solution for problem (|6]) and call it maximum admissible set. 
Notice that the solution for (|6]) might not be unique. The second stage minimizes the total 
transmission power required to support the admitted links: 

, min Y.k(LS, Pk 

S.t. P (SINRfc(p) > 7fc) > 1 - e, ke So, ^^'> 
0<Pk<Pk, k e So- 
Due to the choice of So, power control problem Q is feasible. 

IV. Sample Approximation and Reformulation 

In general, the chance SINR constrained optimization problem ([6]) and problem (|7]) are difficult 
to solve exactly, since it is difficult to obtain the closed-form expression of ([3]) (even though 
CDI is known). In this section, we first approximate the computationally intractable chance 
SINR constraints via sampling, and then reformulate the sample approximation of problems ([6]) 
and (|7]) as a group sparse optimization problem. Three distinctive advantages of the sample 
approximation scheme are as follows. First, it does not require explicit CDI, but only needs 

'Here we do not assume their distribution. 
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samples of channel realizations, which are already available as part of the operation of current 
wireless systems. Second, the sample approximation technique significantly simplifies problems 
(|6]) and (|7]) by replacing the difficult chance SINR constraint with finitely many simple linear 
constraints. Last but not the least, solving the sample approximation problem returns a solution 



to the original chance constrained problem with guaranteed performance [21 1, |22|. 



A. Sample Approximation 



J 1 J 1 f > N 

We handle the chance SINR constraint via sample approximations [21 1, [25 [. Suppose {g^j, Vkf^-i 
are independent samples drawn according to the distribution of {gk,j,Vk}, we use 

alkPk 



SINR^(p) := 



>7fc, neAr:= {1,2,--- ,iV} 



(8) 



9kjPj 



to approximate the chance SINR constraint (|5]). Intuitively, if the sample size N is sufficiently 
large, then the power allocation vector p satisfying ([8]) will satisfy the chance SINR constraint 
Q with high probability. As shown in [ |22| , [ |23| , we have the following Theorem. 
Theorem 4.1: If 




1 + In - + 





2(ir- l)ln- + ln- 



5 



6 



(9) 



for any 5 G (0, 1), then the solution to problem 



mm e p 
p 



s.t. SINR^(p) > 7fc, keJC, neX, 
< p < p, 

will satisfy each of the chance SINR constraint ([5]) with probability at least (1 —5). 



(10) 



Theorem 4.1 shows that it makes sense to call link k is supported if all sampled SINR 
constraints in ([8]) are satisfied. In fact, the above sample approximation technique admits the 
following explanation. Assume that the channel gains {gkj} and the noise power {77^} are random 
variables (but their distribution is unknown). We observe the channel and noise realizations of all 
links at time slots 1, 2, ■ ■ ■ , A^, and design a power allocation strategy (i.e., support some links 
with appropriate power allocation and remove the others) based on the observed channel and 
noise realizations at these time slots. Then this observe-and-design power allocation strategy can 
be used at the following time slots, and it is capable of satisfying the supported links' chance 
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SINR requirements with high probability if the number of observed time slots is large. The 
equation (|9]) reveals how the probability parameter 5 and the number of observed time slots N 
are related to each other. 

We remark that if problem pO]) is feasible, it has the following property. 



Lemma 4.1 (Balancing Lemma): Suppose p is the solution to problem pO] ). Then for each 
k E IC, these must exist an index I < Uk < N such that SINR^*" (p) = 7^. 



It is easy to show Lemma 4.1 by contradiction. Assume that there exists an index k E JC 



such that SINR^ (p) > 7^ for all n E N . Then we can appropriately decrease pk and strictly 



decrease the objective function of problem ( 10) without violating the other links' sampled SINR 



constraints. The Balancing Lemma induces the following distributed Algorithm 1 to solve the 



sampled power control problem ( fTO| ), which is an extension of the Foschini-Miljanic algorithm 
||3|. In fact, when = 1, Algorithm 1 reduces to the Foschini-Miljanic algorithm. The conver- 
gence of Algorithm 1 can be similarly shown as in [j3|. 



Algorithm 1: A Distributed Algorithm for Problem (10) 



Step 1. Initialization: Given a tolerance e > and a feasible point 

p(0) = (pi(0),P2(0), ■ ■ ■ ,PKm^ , and set t = 0. 
Step 2. Iteration: Update the power allocation vector by 

P'^' + = < mm{SINR^(p(t))}^'=(')'^"^ ^ ' ^ ^ 

Step 3. Termination: If ||p(t+l)— p(t) || < e, then stop; else set t = t+1, 
and go to Step 2. 



B. Sampled Channel Normalization 

To facilitate the reformulation of the sample approximation of problems (|6]) and (|7]) and the 
development of efficient algorithms, we normalize the sampled channel parameters. To this aim, 
we use 



Ck 



alkPkdlkPk ' 9k±Pk 




to denote the normalized noise vector of link k. Define 

f 1, if«: = i; 

I 9lkPk 
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and 

A,= (ai;a^;---;af)GM^x^, G /C. 

Notice that the entries of the /c-th column of are one, and all the other entries are nonpositive. 
This special structure of (A; G /C) will play an important role in the reformulation of the 
sample approximation problem. Furthermore, we let 

c = (ci;c2;---;c^)GM^^^^ 

and 

A = (Ai;A2;-- - ; A;,) G M^^^^. 
With the above notations, we can see that 

SINR^(p) > 7fc, V n G AT ^ Afcq > Cfc. 

Consequently, problem (|6]) can be equivalently rewritten as 

max 151 

s.t. Afcq- Cfe > 0, G 5 C /C, (11) 
< q < e, 

and problem Q can be restated as 

min ^ pkqk 

s.t. Afcq - Cfc > 0, A; G 5o, ^^^^ 

< gfc < 1, k e So, 

where So is the solution of the sampled admission control problem ( fTT] ). It is easy to see that the 



sampled power control problem ( 12) is feasible and thus can be efficiently solved by Algorithm 1. 



C. Group Sparse Minimization Reformulation 



It turns out that the two-stage sampled joint power and admission control problem ( [TT] ) and 



( [T2| ) can be reformulated as a single-stage group sparse optimization problem 

min Xlfc II "^ax{cfc - Afeq, 0} llo + ap^q 



q ■ (13) 
s.t. < q < e, 
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where 



< a < «! := 1/e p. 



(14) 



Actually, if there are more than one maximum admissible set (i.e., the solution for problem 
( [TT] ) may not be unique), the formulation ( [T3| ) is capable of picking the one with minimum total 



transmission power as a result of the second term in the objective of ( [13] ). We conclude the 



above results in Theorem |4.2| and relegate its proof to Appendix |Aj 

Theorem 4.2: Suppose q* is the solution to problem ( fT3] ). Then the optimal value of problem 
( fT3| ) lies in (M, M + 1) if and only if the optimal value of problem ( fTT] ) is equal to M (i.e., 
there exists M C K, with \M\ = M such that A^q* > for all k G M). In fact, the set of 
links indexed hy { k E JC \ A^q* > } (whose cardinality is M) is simultaneously supportable 
by the power allocation p* = p o q*. Moreover, p^q* is the minimum total transmission power 
required to support any M links in the network. 



The formulation ( fT3| ) also has the following property stated in Proposition 4.1 which is mainly 
due to the special structure of A^. The proof of Proposition 4.1 can be found in Appendix |B] 



Proposition 4.1: Suppose that q* is the solution to problem ( [13] ) and link k is supported at 
the point q* (i.e., — A^q* < 0). Then there must exist an index 1 < n^^ < N such that 
(Cfe - Afcq*) = 0. 



Proposition 4.1 implies that problem ( [13] ) can be viewed as an (nontrivial) extension of problem 
(]3]). In fact, we know from Proposition 4.1 that when = 1, the solution of problem ( [T3] ) satisfies 
c — Aq* > (i.e., (c — Aq*)^ = for supported links and (c — Aq*)^ > for unsupported 
hnks), and problem ( [T3] ) reduces to problem ([3]). Since problem (]3]) is strongly NP-hard to solve 
to global optimality and strongly NP-hard to approximate to constant factor global optimality 
[|4|' [ fTO| > it follows that problem ( [T3] ) is also strongly NP-hard to solve and approximate. 

A main difference between problem ( [T3] ) and problem ([3]) is as follows. In problem ([3]), if link 
k is supported, then — A^q* in this case is a scalar and equals zero; while in problem ( [T3] ), 
if link k is supported, — A^q* is not necessarily equal to zero. This is the reason why we 
do not just put Cfe — A^q in the objective of ([13]) but put max {c^ — A^q, 0} there. To illustrate 
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A = (Ai;A2) 



/ 1 

1 

-0.3 



-0.5 
1 









/ \ 

1 ^-^ \ 




0.5 


, c = 






0.5 




[0.5) 



It is simple to check that the only possible way to simultaneously support the two links {1, 2} 
is setting q = e. However, 



Aie — Ci 




and A2e — C2 




V. Efficient Deflation Approaches for the Sampled JPAC Problem 
In this section, we develop efficient convex approximation-based deflation algorithms to solve 



the sampled JPAC problem ( [13] ). As we see, problem ( [13] ) has a discontinuous objective function 
due to the first term. However, it allows for an efficient convex relaxation. We first relax problem 



( [13] ) to a convex problem, which is actually equivalent to an SOCP, and then design efficient 
algorithms for solving the relaxation problem. The solution to the relaxation problem can be used 
to check the simultaneous supportability of all links in the network and guide an iterative link 
removal procedure (deflation). We conclude this section with two convex approximation-based 



deflation algorithms for solving the sampled joint control problem ( [13] ) 



A. Convex Relaxation 



Since problem ( [13] ) is strongly NP-hard, we consider its convex relaxation problem 

min ||max{cfc - Ayfcq, 0}||2 + ap^q 
q 

s.t. < q < e. 



To show the convexity of problem ( [T5] ), it suffices to show that ||max {c^ — A^q, 0] 
with respect to q. In fact, for any qi, q2, and A G (0, 1), we have 

||max{Cfc - Ak (Aqi + (1 - A)q2) ,0}||2 
= ||max{A(cfc - A^qi) + (1 - A) (c^ - Afcq2) ,0}||2 

< ||Amax{cfe - A^qijO} + (1 - A)max{cfe - Afcq2,0}||2 

< A ||max{cfc - Afcqi,0}||2 + (1 - A) ||max{cfc - A^qsjOjUg , 



(15) 



IS convex 



(16) 
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where the last inequality comes from the convexity of 



Compared to problem (13), the objective function of problem (15) is continuous in q, but 



is still nonsmooth. We give the subdifferential [26| of the function ||max {c,fc — A^q, 0}||2 in 



Proposition |5.1[ which plays an important role in the following analysis and algorithm design. 
The proof of Proposition 5.1 is provided in Appendix |Cj 

Proposition 5.1: Define /ife(q) = || maxjcfc — A^q, 0} ||2. Then 

5/^fc(q) = {-Ars|s>0,||s||2<l}, 

where — Ajtq = 0. Further, denoting Mk^ = {n \ (c^ — Afcq)„ > 0} 7^ 0, 



V/ifc(q) 



E„6A/-+ (Cfc - Afcq)„ _Al max {c^ - A^q, 0} 



II max {Cfc - Afcq, 0} II2 || max {c^ - A^q, 0} II2 ' 
One may ask why we do not use a smooth convex relaxation instead of ([15])? The reason is 



that, with the use of the nonsmooth mixed £2/^1 norm, the formulation ( [T5| ) is able to characterize 
and induce the group sparsity of the vector x := max {c — Aq, 0} (recall that the objective of 



problem ( 13 ) actually is to find a feasible q such that the vector x is as sparse as possible in the 



group sense). To understand this, observe that the first term of the objective function of ( [15] ), 
the £1 norm of the vector (||xi||2, IIX2II2, ■ ■ ■ , ||x/^||2)"^ , is a good approximation of its £0 norm, 
which is equal to the £0 norm of the vector (||xi||o, ||x2||o, ■ ■ ■ , ||x/c||o)"^ • More discussions on 
using the ^2/^1 norm to recover the group sparsity can be found in [[27l. 



We now discuss the choice of the parameter a in ( [15] ). The parameter a in ( [15] ) should be 
chosen appropriately such that the following "Never-Over-Removal" property is satisfied: the 



solution of problem ( [15] ) should simultaneously support all links at their desired SINR targets 
with minimum total transmission power as long as all links in the network are simultaneously 



supportable. Otherwise, since the solution of ( [T5[ ) will be used to check the simultaneous 
supportability of all links and guide the links' removal, it may mislead us to remove the 



links unnecessarily. Theorem 5.1 gives an upper bound of the parameter a to guarantee the 



"Never-Over-Removal" property. The proof of Theorem 5.1 is mainly based on Lemma 4.1 and 
Proposition 5.1 (see Appendix [D|). 

Theorem 5.1: Suppose there exists some vector q such that < q < e and Aq > c. Then 



the solution of problem (15) with 



< a < ^2 



min {c} 
max {p} 



(17) 



February 26, 2013 



DRAFT 



13 



can simultaneously support all links at their desired SINR targets with minimum total transmis- 
sion power. 



Combining ( fT4| ) and ( [17] ), we propose to choose the parameter a in ( [15] ) according to 

a = min {ciOi, 0202} , (18) 

where Ci,C2 G (0, 1) are two constants. 

Link Removal Strategy. The solution of problem ( [T5] ) can be used to guide the link removal 
process. In particular, by solving ( [T5| ) with a given in ( [TSj ), we know whether all links in the 
network can be simultaneously supported by simply checking if its solution q satisfies Aq > c. 
Furthermore, if all links in the network can not be simultaneously supported, we need to remove 
at least one link from the network. In particular, picking the worst sampled channel index 
fik = argmax{c,fc — A^q} and letting ql = maxjc^ — A^q} for /c G /C, we may remove the 
link with the largest interference plus noise footprint [jTj| 

k = arg max <^ ^ kfej l^j + \^]!k\Qk + C > , (19) 
I j^k j^k J 

the link with the largest excess transmission power [|4| 

k = arg max <^ ^ | a^'^ Ifj + ^l I ^fc p (20) 
I j^k j^k ) 

or the link with the maximum feasibility violation 

/c = argmax{|| max{cA: - Afcq} II2} . (21) 



We remark that the above link removal strategies ( [T9[ ) and ( [20[ ) are different from the ones 
given in ||T|, [|4| due to the presence of channel normalization. We compared the three removal 



strategies ( [19] ), ( [20| ), and (|2Tj) in our simulations, and found that the first one works slightly 
better than the second one (in terms of the number of supported links and the total transmission 
power), and both of them works much better than the last one. Therefore, we shall use the 



removal strategy (19) in the proposed deflation algorithms. In the next subsection, we design 



efficient algorithms to solve the convex but nonsmooth problem (15). 
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B. Optimal Solution for Relaxation Problem ( [T5| ) 

By introducing auxiliary variables x = (xi;x2;-- - ; xx) and t = (ti;t2;--- problem 
( [T5| ) can be transformed into the following SOCP 

min tk + ap'^q 

q,x,t 

s.t. Ilxfcll < tfc, ke IC, 

c — Aq < X, (22) 

< X, 

< q < e, 



which is convex and smooth and can be solved by using the standard solver like CVX |28|. 



Problem ( [15] ) can also be solved directly by the subgradient projection method p9|. The 



subgradient of the objective of problem ([15]) can be easily obtained from Proposition 5.1[ and 



the projection of a point q onto the box Q = {q | < q < e} admits a closed-form solution, 



I.e., 



ProjQ(q) = min {e, max {q, 0}} . 

With an appropriate step size rule, the subgradient projection method is globally convergent. 
However, it has a slow convergence rate that is typically sublinear O (l/v^) , where t is the 
iteration index, and is sensitive to the choice of the initial step size. 



Next, we develop a custom-design algorithm for problem ( [T5[ ) by first smoothing the problem 
and then applying the efficient projected alternate Barzilai-Borwein (PABB) algorithm [30|, [[3T|. 
Specifically, we smooth problem ([T5]) by 



mm /^(q) = X;fc Yl|max{cfe-Afcq,0}||2 + /u2 + ap^q ^^^^ 
s.t. < q < e, 



where /i > is the smoothing parameter. By Proposition 5.1 the objective function /^(q) of 



problem ( [23] ) is differentiable everywhere and its gradient is given by 

f \ -A^maxjcfc - A^q, 0} 
v/^iqj - 2^ , 

keK W||max{cfc - A^q, 0}||2 + /i^ 



As the parameter fi tends to zero, /^(q) uniformly approaches the objective function in ([T5[) and 



hence the solution of the smoothing problem (23) converges to the solution of problem (15). 
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Now we apply the PABB algorithm pO| , pT[ | to solve the smoothing problem ([23]). The 
projected Barzilai-Borwein algorithm is developed to solve the optimization problem 

min /(x) 

X 

s.t. X G X, 

where /(x) is dilferentiable and the projection onto the feasible set X is relatively easy to 
compute. At the t-th iteration, the algorithm updates the design variable x by 

x(t + 1) = Projx [x(t) - aBBWVx/(x(t))] , 

where the step size «BB(i) is either 

with s(t) = x(t) - x(t - 1) and y(t) = Vx/(x(t)) - Vx/(x(t - 1)). In particular, the PABB 
algorithm alternately uses the large step size aBBi and the small step size Q!bb2- 



When using the PABB algorithm to solve problem (23 1, we employ the continuation technique. 



That is, instead of using a tiny /i, we solve ([23]) with a series of values for /i, which are gradually 
decreasing, to obtain an approximate solution of ( [15] ). It turns out the continuation technique 
can reasonably improve the computational efficiency. 

C. Convex Approximation-Based Deflation Algorithms 

The basic idea of the proposed convex approximation-based deflation algorithm for the sampled 



JPAC problem ([13]) is to solve power control problem ( [T5[ ) and check whether all links can be 
supported or not; if not, remove a link from the network, and solve a reduced problem ( ]T5] ) 
again until all the remaining links are supported. 

As in [[10|, to accelerate the deflation procedure (avoid solving too many optimization problems 



in the form of ([T5[)), we derive an easy-to-check necessary condition for all links in the network 
to be simultaneously supported. If all links can be simultaneously served, there exists a vector 
q such that < q < e and Aq > c. By the definition of A, we have 

q > (max {ci} ; max {02} ; • • • ; max {cx}) := c™'^''. 

Define 

^ = A^e 
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and set /x+ = max {/x, 0} and ^_ = max {—fi, 0}. It is obvious that fi = fi^ — /j, . Multiplying 



from both sides of Aq > c, we get that (/x+ — /x_)^ q > e^c. Moreover, we can obtain 

M+e > Ai+q > Ai^q + e^c > /x^c""'''' + e^c, 

where the first inequality is due to q < e and the last one is due to q > c^^^. Therefore, the 
condition 

/x^e - (//^c'"'^" + e^c) > (25) 
is necessary for all links to be simultaneously supported. 



We can use the necessary condition ( [25| ) in the link removal process. In particular, if (|25j) is 
violated, we should drop at least one link from the network. We propose to remove the link ko 
according to 



ko = arg max <^ \ak,j \ + |aj,fc | + L (26) 

k^fC. I I 

L j^k j^k ) 

which corresponds to applying the SMART rule Q to the normalized sampled channel and 



substituting q = e. In ( [26| ), a^j and Ck are the averaged sample channel gain and noise, i.e., 

O'k — afc,2, ■ • • , O'k.K) — , Cfc — ^ , ft fc yV. 



The proposed convex approximation-based deflation algorithm framework for problem ( [13] ) is 
described in Algorithm 2. In the framework, if the power control problem ( fT5| ) is solved via the 
SOCP formulation ( [22] ), we call the corresponding algorithm SOCP-D; while if problem ( |T5| ) 



is solved via the smoothing problem ([23]) and the PABB algorithm, we call the corresponding 
algorithm PABB-D. The SOCP-D algorithm is of polynomial time complexity, i.e., it has a 
complexity of 0{N'^■^K^^^{N + K)), since it needs to solve K SOCP problems (|22|) and solving 



one SOCP problem in the form of ([22|) requires 0{N^-^K'^-^{N + K)) operations [|32| Page 423]. 
We should remark that the postprocessing step (Step 5) in the framework aims at admitting the 
hnks removed in the preprocessing and admission control steps. Specifically, we enumerate 
all the removed links and admit one of them if it can be supported simultaneously with the 
already supported links. If there are more than one such candidates, we pick the one such 
that the minimum total transmission power is needed to simultaneously support it with the 
already supported links. The postprocessing step is terminated if no such candidate exists. The 
specification of the postprocessing step can be found in Appendix |Ej 
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Algorithm 2: A Convex Approximation-Based 
Deflation Algorithm Framework 
Step 1. Initialization: Input data (A,c,p) . 



Step 2. Preprocessing: Remove link Ajq iteratively according to ([26]) 
until condition (1251) holds true. 



Step 3. Power control: Compute parameter a by ( 18 1 and solve problem 



( [15] ); check whether all links are supported: if yes, go to Step 5; else 
go to Step 4. 



Step 4. Admission control: Remove link according to ( ]T9] ), set 
K, = K,/ {ko} , and go to Step 3. 

Step 5. Postprocessing: Check the removed links for possible admis- 



sion. 



VI. Numerical Simulations 

To illustrate the effectiveness and efficiency of the two proposed convex approximation-based 
deflation algorithms (SOCP-D and PABB-D), we present some numerical simulation results in 
this section. The number of supported links, the total transmission power, and the execution 
CPU time are used as the metrics for comparing different algorithms. 

Simulation Setup: First, we generate the nominal channel parameters as in [[T|, i.e., each 
transmitter's location obeys the uniform distribution over a 2 Km x 2 Km square and the location 
of each receiver is uniformly generated in a disc with center at its corresponding transmitter and 
radius 400 m, excluding a radius of 10 m. The nominal channel gains are given by 



9k, 



where d^j is the Euclidean distance from the link of transmitter j to the link of receiver k. The 
nominal noise power is set to be 



Vk 



-90 dBm, W keJC. 



Second, we generate samples of the channel parameters g^j and 77^ for n 
Specifically, for n = 1, we set 



1,2,- 



.N. 



9lj = 9k,j, Vk = Vk, k,j e /C; 



(27) 
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for n = 2, 3, ■ ■ ■ , iV, we set 




(28) 

7^- = fi,*{l + S- 



In the above, and (J^ are samples of the corresponding random variables ^kj and (k, both 
of which obey the standard Gaussian distribution, i.e., ^kj ~ A/'(0, 1) and (^k ~ A/'(0, 1); ^'"^j^ = 
max„=2,3,- ,Af {l^fc.jl} ^iid Cr'''' = max„=2,3,- ,Af {ICfc 1} ; and S G [0, 1) is a given number, which 
guarantees that all generated g^j and r^^ are strictly greater than zero. Notice that the parameter 
S controls the magnitude of the perturbation added to the nominal channel parameters. The 
larger the parameter S is, the larger the perturbation is. 

In addition, each link's SINR target is set to be 7,^ = 2 dB (y k E /C) and the power budget 
of the link of transmitter k is set to be = 2p^ (V /c G /C), where is the minimum power 
needed by link k to meet its SINR requirement without any interference from other links under 
the nominal channel condition. 

Benchmark: From the way of generating channel parameters, we know that if g^j = gkj and 
Vk — — 1) 2, ■ ■ ■ , A^, then the number of supported links in this case (without 

perturbation) should be greater than or equal to the number of supported links under the channel 



conditions with perturbation given in p7] ) and ( [28] ). In addition, if the number of supported 
hnks under these two cases are equal to each other, the total transmission power in the former 
channel condition should be less than the one in the latter channel condition. In particular, when 



aU samples of the channel parameters are equal, the corresponding sampled JR^C problem (13) 



reduces to problem ([3]), which can be solved efficiently by the NLPD algorithm in [jTOj. The 



solution given by the NLPD algorithm will be used as the benchmark to compare with the two 
proposed algorithms, since the NLPD algorithm was reported to have the close-to-global-optimal 



performance in [10|. For fair comparison, we choose the parameter a according to ([18]) and use 



the removal strategy ( [T9[ ) in the NLPD algorithm as in the two proposed algorithms. 

Choice of Parameters: We set the parameters e, 5, and i^' in ([9]) to be 0.1, 0.05, and 10, 
respectively. We remark that K in equation ([9]) is the number of supported links but not the 

^This corresponds to the special case where S — 0. 
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TABLE I 

Statistics of the Number of Supported Links of 200 Monte-Carlo Runs. 



Number of Total Links 


Parameter 5* 


Algorithm 


Statistics of the Number of Supported Links 


4 





Benchmark 


664=2*20-^3*96-1-4*84 


4 


0.1 


SOCP-D/PABB-D 


647=2*26-1-3*101-1-4*73 


4 


0.2 


SOCP-D/PABB-D 


619=1*2-1-2*34-1-3*107-14*57 


12 





Benchmark 


1468=5*8-1-6*35-1-7*65-^8*68-^9*21-1-10*3 


12 


0.1 


SOCP-D/PABB-D 


1398=4* 1-1-5* 13-1-6*49-1-7*75-1-8*5 1-h9*8-h10*3 


12 


0.2 


SOCP-D/PABB-D 


1291=4*5-1-5*27-1-6*78-1-7*57-1-8*28-^9*5 


20 





Benchmark 


1988=6*l-H7*5-i-8*21-H9*49-t-10*61-i-ll*39-Fl2*16-Hl3*5-Hl4*3 


20 


0.1 


SOCP-D/PABB-D 


1841=5*1-1-6*3-1-7*17-1-8*39-^9*63-^10*40-1-1 1*27-Fl2*7-I-13*3 


20 


0.2 


SOCP-D/PABB-D 


1676=6*8-1-7*44-^8*66-^9*43-1-10*27-^11*9-^12*3 


28 





Benchmark 


2345=9*11-1-10*25-^11*61-1-12*48-1-13*31-1-14*19-1-15*2-1-16*1-1-17*2 


28 


0.1 


SOCP-D/PABB-D 


2189=8*4-1-9*23-^10*54-1-1 l*56-Hl2*39-i-13*13-Hl4*8-i-15*3 


28 


0.2 


SOCP-D/PABB-D 


1975=7*8-1-8*23-1-9*51-1-10*56-^1 1*34-1-12*23-1-13*4-1-14*1 



number of total links. Substituting these parameters in (|9]), we obtain N* = 200, and we set 



= 200 in all of our simulations. Both of the parameters ci and C2 in ( [18] ) are set to be 0.999. 



When applying the PABB algorithm with the continuation technique to solve the power control 



problem ( [15] ), we always set the smoothing parameter /x to be 1 at first, and then gradually 
decrease it by multiplying a factor of 0.1 until it becomes less than 10^^. Two values are used 



for the S in ([28]), which are 0.1 and 0.2. Finally, we use CVX [[28J to solve the SOCP problems 
in the SOCP-D algorithm. 

Simulation Results and Analysis: Table [ij summarizes the statistics of the number of supported 
links of 200 Monte-Carlo runs in two cases where 5 = 0.1 and 0.2, respectively. For instance, 
"664 = 2* 20 + 3* 96 + 4* 84" in the fourth column of Table | stands for that when K = 4, 664 
hnks are supported in these 200 Monte-Carlo runs, and amongest them, 2 links are supported 
20 times, 3 links are supported 96 times, and 4 links are supported 88 times. Figs. [TJ [2] and 
[3] are obtained by averaging over the 200 Monte-Carlo runs. They plot the average number of 
supported links, the average total transmission power, and the average execution CPU time of 
the proposed SOCP-D and PABB-D algorithms and the benchmark versus different number of 
total links, respectively. 

We can observe from Fig. [T]that the number of supported links by the two proposed algorithms 

Febmary 26, 2013 DRAFT 



20 




12 16 20 

Number of Total Links 



Fig. 1. Average number of supported links versus tiie number of total links with 5 = 0.1 and 0.2. 




12 16 20 

Number of Total Links 



Fig. 2. Average total transmission power versus the number of total links with S* = 0.1 and 0.2. 



is slightly less than the benchmark, which serves as an upper bound of the number of supported 
links by the proposed algorithms. Table |I} Fig. [1] and Fig. [2] show that the two proposed 
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Fig. 3. Average execution time (in seconds) versus the number of total links with S — 0.1 and 0.2. 



algorithms always return the same solution to the sampled JPAC problem ( [13] ), i.e., supporting 
same number of links with same total transmission power. However, Fig. [3] shows that the 
PABB-D algorithm substantially outperforms the SOCP-D algorithm in terms of the average 



CPU time. This is not surprising, since in each power control step (i.e., solving problem ([15])), 
the custom-design algorithm is used to carry out power control in the PABB-D algorithm while 
a general-purpose solver CVX is used to update power in the SOCP-D algorithm. 

We know the fact that when we add the perturbation to the nominal channel parameters (cf. ([27]) 



and ([28])), either the number of supported links under the perturbed channel condition equals the 
number of supported links under the nominal channel condition but the total power transmitted 
in the former case is greater than the one in the later case, or the number of supported links in the 
perturbed channel condition is less than the one in the nominal channel condition. As plotted in 
Fig. [2] as the number of total links in the network increases from 24 to 28, the difference of the 
average total transmission power by the benchmark and the two proposed algorithms increases. 
This, combining the above fact, implies the difference of the number of supported links by the 
benchmark and the two proposed algorithms will become large as the number of links in the 
network increases from 24 to 28. This is consistent with the simulation result in Fig. [T] 
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Comparing the two scenarios where 5 = 0.1 and S = 0.2, it can be seen from Fig. [T]that the 
number of supported links by the two proposed algorithms in the scenario 5 = 0.2 is less than 
the one in the scenario S = 0.1. This makes sense because the magnitude of the perturbation 
in the scenario S* = 0.1 is generally smaller than the one in the scenario S = 0.2. The fact that 
the number of supported links in the scenario S = 0.1 is larger than the one in the scenario 
S = 0.2 also reveals the reason why the average total transmission power in the former case 
is larger than the one in the latter case as shown in Fig. [2j Notice that in Fig. [2[ the average 
total transmission power in the scenario S* = 0.1 is smaller than the one in the scenario S* = 0.2 
when K = A. This is because in this case the average number of supported links in the two 
scenarios are nearly equal to each other. We also point out that the execution CPU time of 
the two proposed deflation algorithms mainly depends on how many times the power control 



problem ( 15 1 is solved. In general, the larger number of links are supported, the smaller number 



of links are removed from the network and the smaller number of power control problems in 



form of (15) are solved. Therefore, the average CPU time of the proposed algorithms in the 



scenario S* = 0.1 is larger than the one in the scenario S* = 0.2; see Fig. [3] 

VII. Conclusions 

In this paper, we have considered the chance SINR constrained JPAC problem, and have 
proposed two sample approximation-based deflation approaches for the problem. We first ap- 
proximated the computationally intractable chance SINR constraints by sampling, and then 
reformulated the sampled JPAC problem as a one-stage group sparse £o-iTiinimization problem. 
Furthermore, we relaxed the NP-hard ^o-minimization problem to a convex i*! -minimization 
problem (equivalent to an SOCP) and used its solution to check the simultaneous supportability 
of all links and guide an iterative link removal procedure (deflation), resulting in efficient deflation 
algorithms (SOCP-D and PABB-D). 

The proposed approaches are particularly attractive for practical implementations for the 
following reasons. First, since the proposed approaches are based on sample approximations, 
they do not require explicit CDI but only need samples of channel realizations. Consequently, 
the proposed approaches are distributionally robust. Second, the two proposed approaches enjoy 
a low computational complexity. The SOCP-D approach has a polynomial time complexity. 
To further improve the computational efficiency, the special structure of the SOCP relaxation 
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problem is exploited, and an efficient algorithm, PABB-D, is custom designed for solving it. The 
PABB-D algorithm significantly outperforms the SOCP-D algorithm in terms of the execution 
CPU time. Finally, our simulation results show that the proposed approaches are very effective 
by using the nonrobust NLPD algorithm as the benchmark. 
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Appendix A 



Proof of Theorem 4.2- It is simple to see the equivalence between problem ( [TT] ) and the 
intermediate problem 

min Xlfc II ~ Afe^'O} llo 



1 ^" (29) 
S t. < q < e 



in the sense that the optimal value of problem ([TT]) is M if and only if the minimum value of 



problem (29) is - M. 



We now establish the equivalence between problems ([T3]) and (29) (in the sense of supporting 



the maximum number of links). We claim that the optimal value of ( [29| ) is M if and only if the 



optimal value of ([T3]) is in the interval (M, M + 1). We argue the "only if and "if directions 



separately. Recall the fact < a < ai = l/p"^e (cf. ([14])) which implies 

a p^q < 1, for any < q < e. 

Consequently, the total contribution from the second term in the objective function of ( [T3] ) cannot 
exceed 1, regardless of the power allocation q. This immediately shows that if the optimal value 
of ( [29] ) is Af, then the optimal value of ( [T3] ) must be in the interval (M, M + 1). To argue 



the converse, we note that the function Xl/tG/c II Hiaxjcj. — A^q, 0} ||o is discontinuous with an 
increment of 1, implying that 

II max {Cfc - Afcq, 0} ||o + a p^q 

k&K 

> J^||max{cfc-Afcq,0}||o (30) 

> M+l>M + ap^q, 
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for any choice of < q < e, as long as J2keK. II ^^-^ {^^^ ~ A^q, 0} ||o > M. Thus, the global 



minimum of ( [13] ) must be achieved at a power allocation for which 

^ II max {Cfe - Afcq, 0} ||o = M 

keK. 

holds, i.e., J2keJC II ^^^^ ~ A^q, 0} ||o is fully minimized by ([29]). This establishes the equiv- 
alence between ( [29] ) and ( [T3] ). 



Finally, if there are multiple sets of M links that are simultaneously supportable, then they all 
induce the same objective value in the first term of the objective function in ( [T3] ). In this case, 
the second term (i.e., p-'^q) will play the role to select the one set of M links which requires the 
least amount of total transmission power. In other words, p^q* = e^p* is the minimum total 
power required to support any M links in the network. 

Appendix B 



Proof of Proposition 4.1- We prove Proposition 4.1 by contraction. Assume that link k is 



supported in ( [13] ) but A^q* > holds true. Then we can construct a feasible point q satisfying 
II max {Cfc - Afeq, 0} ||o + ap'^q < ^ || max {c^ - A^q*, 0} ||o + ap^q*. (31) 

k k 

Define q = (gi, g2, ■ " " , Qk)^ with 

max {(Efc - Afc)q* + c^} , if j = k; 
q*, if J ^ k. 

Recalling the definitions of and A^, we know — A^ is a nonnegative matrix, and thus 

Qk = max {(Efc - Afc)q* + Cfc} > 0. 

Since 

Afcq* = Efcq* - (Efc - Afc)q* = g*e - (E^ - Afe)q* > c^, 

it follows that 

ql > max {(Efc - Afc)q* + c^} = qk- 

Hence, q is feasible (i.e., < q < q* < e and the inequality q < q* holds true strictly for the 
A;-th entry) and 

p'^q < p^q*. (32) 
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Moreover, it follows from the definition of that A^q, > c^. For any j ^ k, if A^q* > Cj, 
we have Ajq — Cj > A^q* — Cj > 0, where the first inequality is due to the fact that all entries 
of Aj except its j-th column are nonpositive and the fact q < q*. Consequently, there holds 
J* C J, where J* = {j \ Ajq* > Cj} and J = {j \ A^q > Cj} . Thus, we have 

Y^kejc II "^^^ {^k - Afcq*, 0} llo 
= II Kiax{cfe - Afcq*,0} llo 

(33) 

> Efc^j II inax{cfe - Afcq,0} llo 

= Y.k&K.W'^^^i^k- AkCi,0}\\o. 
Combining ( [32] ) and ( [33] ) yields pT] ), which contradicts the optimality of q*. This completes the 



proof of Proposition 4. 1 



Appendix C 

To prove Proposition |5.1[ we first consider the simple case where h{q) = \\ max{q, 0} ||2. 
Lemma C.l: Suppose /i(q) = || max {q, 0} II2. Then dh(0) = {s | s > 0, ||s||2 < 1} . If there 
exists i such that (q)j > 0, then h(q) is differentiable and 

VMq) = „ '""''''■°'„ . (34) 

II maxjq, 0} II2 

Proof: Define 5 = {s | s > 0, ||s||2 < 1} . We claim dh{0) = S. On one hand, taking any 
s G 5, we have that 

h{q) = II max {q, 0} II2 > max {q, 0} > s^q = h{0) + s^(q - 0), V q, 

where the first inequality is due to the Cauchy-Schwarz inequality and the fact ||s||2 < 1, and 
the second inequality is due to s > 0. This shows that S C dh(0) according to the definition of 
dh(0) p6| . On the other hand, to show dh{0) C iS, it suffices to show that any point s ^ 5 is 
not a subgradient of h{q) at point 0. In particular, if ||s||2 > 1, then 

/i(q) = II max {q, 0} II2 < ||q||2 = 1 < ||s||2 = s^q = /i(0) + s^(q - 0) 

with q = s/||s||2. Thus, the subgradient s of h{q) at point must satisfy ||s||2 < 1. If ||s||2 < 1 
but (s)i < (without loss of generality), we test the point q = (—1, 0, ■ • ■ , 0)"^ , and obtain 

/i(q) = ||max{q,0} ||2 = < - (s)^ = s^q = /i(0) + s^(q - 0). 
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Consequently, the subgradient s of h{q) at point must satisfy ||s||2 < 1 and s > 0. Hence, 
dh{0) = {s|s > 0, ||s||2 < 1}. 

Next, we show that h{q) is differentiable at the point q which has at least one positive entry, 
and the corresponding gradient is given in ([34]). In fact, although the function max{g,0} is 
nondifferentiable at point q = 0, its square f{q) = (max{g, 0})^ is differentiable everywhere; 
i.e., 

' 0, if g < 0; 
2q, if X > 0. 

According to the composite rule of differentiation, we know that the gradient of h(q) is given 
by ( [34] ). This completes the proof of Lemma C.l[ ■ 



Equipping with Lemma C.l , we now can prove Proposition 5.1 Assume that c^. — A^q = 0. 
Then for any q and any s satisfying ||s||2 < 1 and s > 0, we have 

/ifc(q) = II max{cA: - Afcq, 0} II2 

> (cfc — Afcq) (from Lemma ICT]) 

= s'^ (Afcq - Afcq) 

(35) 

(q-q) 

= II max {Ck - Afcq, 0} II2 + (-A^s)'^ (q - q) 
= hk{^) + {-Alsf {q-q) , 
which shows that all vectors in {— A|^s | s > 0, ||s||2 < l} are subgradients of hk{q) at point 



q. In the same way as in the proof of Lemma C.l we can show that if s does not satisfy 



II s II 2 < 1 and s > 0, the inequality in ([35]) will violate for some special choices of q. Hence, 
dhkiq) = {-A^s|s > 0, ||s||2 < 1}. 

If J\fft 7^ 0, we know from the composite rule of differentiation and Lemma C.l that /ifc(q) 
is differentiable and its gradient is given by 



V/ifc(q) 



EngA/-+ (Cfc - Afcq)n «)' 

II max{cfc - Afcq, 0} II2 



-A^maxjcfc - Afcq, 0} 
I max{cfc - Afcq, 0} II2 ' 



This completes the proof of Proposition 5A\ 
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Appendix D 



Proof of Theorem 5.1- Suppose all links in the network can be simultaneously supported (i.e., 
there exists < q < e satisfying Aq > c) and q is the solution to problem 



mm p q 

q 



s.t. Aq - c > 0, 
< q < e. 



(36) 



To prove Theorem 5.1, it suffices to show that q is also the solution to problem (15) with 



a G [0, a2\. Moreover, to show q is the solution to problem ( [T5] ), we only need to show that the 



subdifferential of the objective function of problem ( [15] ) at point q contains [26|. Next, we 
claim the latter is true. 



We first characterize the subdifferential of the objective function of problem ( [T5| ) at point q. 
It follows from Lemma 4.1 that there exists X = {ni,n2, ■ ■ ■ such that q is the solution 



to the following linear system 

Axq:= [ar;a^^--- ;a^-]q 



1 



) ■■= ex. 



Recalling the definition of a)! (see Subsection IV-B), we know that I — Aj is a nonnegative 
matrix. Moreover, from |i33| Theorem 1.15], Aj is nonsingular, Aj^ is nonnegative, and 

< e^Aj^Cx = e^q < K. (37) 



Define /ifc(q) = ||max{cfc — A^q, 0}||2 for A; G /C. It follows from [26, Theorem 23.8] that the 



subdifferential of the objective function of problem (15) at point q is given by 



<^gk + (yp\gk e dhk{q), k e }c \ . 

[keK. J 
According to Proposition 5.1 dhk{q) contain^aW vectors in 

Sk = [sk {al'-f |0 < Sfc < l}. 



(38) 



Therefore, all vectors in 



S = {-A^s + ap I < s < e} = i - ^ Sfe «*)^ + ap | < < 1, A; G /C I 

I keJC ) 



^The subdifferential of /ife(q) at point q is not necessarily equal to Sk in ([38j. This is because that some other entries (except 
the nfe-th entry) of the vector Ck. — A^q might also be zero. 
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are subgradients of the objective function of problem ([T5|) at point q. 



If G 5, the subdifferential of the objective function of problem (15) at point q contains 



1 26 1, which completes the proof of Theorem 5.1 In fact, consider the vector s = «Aj p. It is 
a nonnegative vector (since A^^ is nonnegative), and each of its entries is less than or equal to 



1 as long as a < • This is because 



T 

e s 



< a max {p} e-^Aj^e 
max{p} ^ 
mmjcxj 
max {p} 



< a 

(a) 

< a 

< 1. 



(39) 



min {c} 



K 



where (a) is due to ([37]) and the fact min {c} < min {cj}. Substituting s = aA^^p into S, we 
obtain 

— AjS + ap = — Aj (aAj^p) + ap = 0. 



Thus, eS. 
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Appendix E 



Algorithm 3: Specification of Postprocessing 
Step 1. Initialization: Input the set of supported links S and the set of 

removed links TZ. 
Step 2. Iteration: 
For j e TZ, 
- Solve the linear program 

min p"^q 
q 



s.t. Afcq-Cfc>0, keS[j{j}, (40) 
< q < e, 



and obtain its optimal value Vj . If problem ( [40] ) is infeasible, 
Vj is set to be +oo. 

- If Vj = +00, set n = n/ {j} . 

End (For) 

Step 3. Termination: If 7^ = 0, then stop; else choose 

j* = argmin {vj} , 
set 5 = 5 U {j*} and n = n/ {j*} , and go to Step 2. 
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